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ABSTRACT 

A  rotating  disc  was  designed  to  fail  due  to  time  dependent  mechanisms. 
Several  discs  were  fabricated  and  an  existing  test  rig  was  developed  to  test 
the  discs  at  design  conditions.  Ten  discs  were  tested  at  steady  state 
under  the  design  conditions  of  2300°F  rim  temperature  and  50,000  rpm  for 
periods  of  0,20  hours  to  25  hours.  An  experimental  failure  distribution 
was  obtained  for  the  ten  discs  and  presented  as  reliability  versus  time. 
Three  different  data  bases  (double  torsion,  stress  rate,  and  stress  rupture) 
were  used  to  calculate  reliability  versus  time  and  the  results  were  com¬ 
pared  to  the  experimental  results.  It  was  observed  that  the  best  correla¬ 
tion  with  the  experimental  results  was  given  by  the  stress  rupture  data 
base .  *7' 
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INTRODUCTION 


Since  the  early  sixties  there  has  been  considerable  interest  i^ 
silicon  nitride  as  an  engineering  material.  Sage  and  Histed 
reviewed  some  of  the  applications  for  silicon  nitride  and  suggested 
the  use  of  silicon  nitride  in  gas  turbines.  A  potential  increase  of 
100®C  in  turbine^  stator  blade  service  temperature  was  forecast. 
Glenny  and  Taylor  reviewed  the  material  properties  of  silicon  nitri¬ 
des  and  noted  the  high  strength^  at  elevated  temperatures  of  silicon 
nitride.  Evans  and  Wlederhorn  observed  the  phenomenon  of  static 
fatigue  in  hot  pressed  silicon  nitride,  and  noted  that  this  presented 
3  design  problem  in  the  application  of  the  material  in  highly  stressed 
applications  at  elevated  temperatures. 

The  phenomenon  of  static  fatigue  is  well  known  in  jaetals  where  it 
may  be  caused  by  slow  crack  growth.  Stverns  and  Price  observed  that 
the  crack  growth  rate  in  steel  at  elevated  temperatures  was  correlated 
with  the  slow  crack  growth  power  equation 

V  =  AI^  (1) 


Evans  and  Johnson  suggested  that  this  equation  could  be  applied  to 
ceramics  and  that  the  cause  of  slow  crack  growth  was  grain  boundary 
sliding.  This  is  due  to  the  glassy  phases  in  the  grain  boundary  sof¬ 
tening  at  elevated  temperature  and  allowing  the  crystalline  grains  of 
silicon  nitride  to  slide  relative  to  one  another  at  the  tip  of  the 
crack  where  the  stresses  are  the  highest.  The  sliding  of  the  grains 
causes  the  crack  co  grow  and  continued  aj^plication  of  stress  and  tem¬ 
perature  could  cause  catastrophic  failure  , 

7 

Williams  and  Evans  proposed  the  dcxuble  torsion  test  to  measure 
the  values  of  the  crack  velocity  exponent,  and  the  premultiplier  in 
Eqn.  1.  The  double  torsion  test  is  a  simple  test  in  concept,  that 
loads  a  pre-cracked  specimen  in  a  manner  that  the  stress  intensity 
factor  remains  constant  as  the  crack  grows.  This  simplifies  the  ana¬ 
lysis  and  makes  it  possible  to  determine  A  and  n  from  the  load  versus 
time  curve. 

rt 

Charles  proposed  stress  rate  testing  to  determine  the  crack 
velocity  exponent  and  demonstrated  its  use  on  glass.  In  stress  rate 
testing  MOR  bars  are  fractured  at  two  different  stress  rates  to  obtain 
two  pairs  of  data  -  the  fast  fracture  stress  and  its  corresponding 
stress  rate.  From  these  sets  of  data  the  slow  crack  growth  parameters 
are  calculated. 

Several  experimenters  ha^e  conducted  stress  rupture  tests  on 
ceramic  materials.  Trantina  conducted  stress  rupture  tests  on 
Norton's  NC-132  material,  a  hot  pressed  silicon  nitride.  His  tests 
were  conducted  in  three-point  bending  on  bars^^f  0.9  inch  span  with  a 
cross  section  of  0.1  inch  x  O.l  inch.  Quinn  conducted  stress  rup¬ 
ture  tests  on  several  different  structural  ceramic  materials.  Hit- 
tests  were  conducted  in  four-point  bending  on  a  fixture  with  a  1.5 
inch  outer  span,  a  0.75  inch  inner  span,  and  a  bar  with  a  0.110  inch 
base  .1  1  a  0.085  inch  height. 


Previous  experimenters  conducted  their  tests  on  small  parts  with 
very  small  volumes  compared  to  the  volume  of  typical  gas  turbine 
components.  Previous  tests  were  also  conducted  at  uniform  tem¬ 
peratures,  but  gas  turbine  components  have  large  temperature  gra¬ 
dients,  This  suggested  the  need  for  generating  a  time  to  failure 
data  base  on  a  component  the  size  of  a  typical  gas  turbine  component 
with  temperature  and  stress  distributions  similar  to  gas  turbine  com¬ 
ponents.  This  data  base  could  then  be  used  as  a  test  case  for  eva¬ 
luation  of  I  adels  of  failure  and  computer  programs  used  to  predict 
time  to  failure. 


Figure  1  -  Hot-spin  test  rig 

The  hot  spin  rig  shown  in  Fig,  1  was  selected  for  testing  of  a 
gas  turbine  component  because  it  had  accumulated  many  hours  of  testing 
and  had  shown  that  it  could  test  turbine  di^^cs  at  speeds  and  tem¬ 
peratures  simulating  actual  engine  conditions^  .  The  rig  development 
required  to  meet  the  test  ^c^nditlons  and  the  initial  calibration  test 
were  previously  described  ,  This  report  reviews  the  most  recent 
results  with  emphasis  on  the  determination  of  experimental  reliability 
versus  time  and  the  comparison  to  calculated  reliability  versus  time. 


II.  DESIGN  OF  THE  DISC 


A  finite  element  model  of  the  rotating  components  was  set  up  for 
calculating  the  temperatures  of  the  disc.  The  model  is  shown  in  Fig. 
2.  The  physical  parts  represented  by  the  finite  elements  are  shaded 
and  labeled  for  identification.  The  finite  elements  that  compose  the 
model  are  rectangular  or  triangular  in  shape.  The  material  thermal 
properties  used  in  the  analysis  are  given  in  Table  I. 
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Figure  2  -  Rotating  components  heat  transfer  model 


The  temperatures  within  the  rotating  components  were  determined 
by  the  boundary  conditions  and  the  thermal  conductivities  of  parts  in 
the  assembly.  Heat  was  applied  to  the  front  face  and  rim  of  the  spin 
disc  with  propane  gas  burners.  The  temperature  of  the  front  face  was 
measured  with  a  radiation  pyrometer  and  the  measured  temperatures  were 
used  as  input  to  the  analysis.  The  tie  bolt  was  cooled  with  air  as 
shown  on  Fig.  2.  The  cooling  air  mass  flow  rate  and  associated  heat 
transfer  coefficients  formed  the  boundry  conditions.  A  portion  of  the 
tie  bolt  cooling  air  was  directed  to  the  metal  curvic  adaptor  for 
cooling  and  the  mass  flow  rate  and  heat  transfer  coefficients  formed 
the  boundary  conditions  for  the  exterior  of  the  tie  bolt,  nut,  and 
metal  curvic  adaptor.  The  dowristream  side  of  the  spin  disc  and  the 
exterior  of  the  ceramic  curvic  spacer  were  heated  with  hot  gases  and 
air  swirling  In  the  rig,  thus  forming  another  boundary  condition.  The 
labyrinth  and  bearing  region  on  the  shaft  form  the  last  boundary  con¬ 
dition.  As  part  of  the  model,  radiation  heat  transfer  was  included 
between  the  spin  disc  bore,  ceramic  curvic  spacer  bore  and  the  tie 
bolt. 
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Calculated  temperature  distributions  were  compared  to  results 
from  temperature  sensfilve  paints  and  adjustments  to  the  model  were 
made  to  obtain  agreement.  Temperature  measurements  of  the  air  in  the 
cooling  streams  were  also  used  as  Input  to  the  model  and  as  a  check  on 
the  validity  of  the  model;  for  example  the  temperature  of  the  tie  bolt 
cooling  air  was  measured  as  it  entered  the  labyrinth  and  as  it  left  the 
tie  bolt.  The  temperature  of  the  air  in  the  labyrinth  region  was  used 
as  input  to  the  model  and  the  air  exit  temperature  was  compared  to 
the  calculated  value  as  a  check  on  the  model. 

Figure  3  illustrates  the  finite  element  grid  used  to  calculate 
the  stresses  in  the  spin  disc.  The  temperature  distribution  was 
derived  from  the  heat  transfer  model  and  transferred  to  the  stress 
finite  element  grid.  The  temperature  contour  plot  at  the  test  design 
conditions  is  illustrated  in  Fig.  4.  The  combined  thermal  and  centri¬ 
fugal  stresses,  at  the  test  condition,  are  Illustrated  in  Figs.  5 
through  9.  These  stress  plots  show  the  maximum  principal  tensile 
stress,  the  tangential  stress,  the  radial  stress,  the  axial  stress, 
and  the  maximum  shear  stress. 


N - .67" 


Figure  3  -  Disc  finite  element  grid 
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Figure  4  -  Isotherms  [®F]  in  disc  at  50,000  rpm. 


Figure  5  -  Maximum  principal  stresses  [psi]  in  disc  combining 
thermal  and  centrifugal  stresses  at  50,000  rpm. 


Figure  8  -  Axial  stresses  (psi]  in  disc  combining  thermal 
and  centrifugal  stresses  at  50,000  rpm. 


Figure  9  -  Maximum  shear  stresses  [psl]  in  disc  combining 
thermal  and  centrifugal  stresses  at  50,000  rpm. 


The  previous  figures  Illustrate  the  combined  thermal  and  mechani¬ 
cal  stress  conditions  at  which  the  stress  rupture  tests  were  con¬ 
ducted.  The  mechanical  stresses  due  to  centrifugal  loads  at  50,000 
rpra  were  calculated  with  the  disk  at  a  constant  temperature  of  78*F, 
and  the  stresses  are  plotted  in  Figs.  10  through  14,  These  stress 
plots  show  the  maximum  principal  tensile  stress,  the  tangential 
stress,  the  radial  stress,  the  axial  stress,  and  the  maximum  shear 
stress,  and  may  be  compared  to  the  combined  stresses  of  Figs.  5 
through  9.  The  elastic  properties  used  in  the  stress  calculations  are 
given  in  Table  II. 


Figure  10  -  Maximum  principal  stresses  [psl]  in  disc, 
centrifugal  stresses  at  50,0000  rpm. 


Figure  11  -  Tangential  stresses  [psl]  In  disc,  centri¬ 
fugal  stresses  of  50,000  rpm. 
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Figure  lA  -  Maximum  shear  stresses  (psl]  In  disc,  centri¬ 
fugal  stresses  at  50,000  rpm. 

Three  data  bases:  double  torsion,  stress  rate,  and  stress  rup¬ 
ture,  were  available  for  the  calculation  of  time  dependent  reliabili¬ 
ties.  l^e  first  two  data  bases  are  used  with  a  formula  derived  by 
Paluszny  and  repeated  here  as  Equation  2.  Definition  of  the  symbols 
Is  given  in  Appendix  C. 


(2) 


R  •  exp  - 


n-2 

n 

ffl 

(In-i— )  ”  + 

a  t 

n-2 

®Ef 

Og  B 

During  this  work,  Eqn,  2  was  reduced  to  a  simpler  form  shown  a^^^qn. 
3,  The  detailed  algebraic  steps  are  given  in  the  Interim  report^  . 


m 

n-2 


(3) 


Equation  3  Is  much  simpler  for  computation  and  was  programmed  for 
digital  computer  use  to  calculate  time-dependent  reliabilities  for 
steady  state  temperature  and  loading  conditions.  The  equation  was 
evaluated  for  each  element  of  the  finite  element  grid  of  Fig.  3  using 
the  maximum  principal  stress,  the  Weibull  modulus,  the  crack  velocity 
exponent,  the  value  of  B,  and  the  fast  tractiire  reliability  of  the 
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element.  The  Welbull  modulus,  crack  velocity  exponent,  and  Che  value 
of  B  were  linearly  interpolated  from  a  table  of  these  parameters  ver¬ 
sus  temperatures,  based  on  the  centroldal  temperature  of  Che  element. 

Equation  3  shows  that  when  time  equals  zero  the  time  dependent 
reliability  is  simply  the  fast  fracture  reliability.  The  fast  frac¬ 
ture  reliability  was  evaluated  by  integrating  the  normal  stress  around 
the  unit  sphere  in  each  finite  ^l^m^t^^and  evaluating  the  Weibull 
volume  integral  for  each  element^  ,  ,  .  The  double  torsion  and 
stress  rate  data  bases  yield  Che  parameters  "B"  and  “n"  essential  for 
the  evaluation  of  Eqn.  3.  Figure  15  shows  the  calculation  flow  path 
utilizing  double  torsion  data.  Fast  fracture  testing  produced  a 
collection  of  fast  fractu^l  data  which  was  analyzed  with  maximum  like¬ 
lihood  estimation  methods  to  determine  the  characteristic  modulus  of 
rupture  and  the  Welbull  modulus.  The  stress  and  temperature  distribu¬ 
tions  from  finite  element  analysis  along  with  the  Weibull  parameters 
were  input  Co  a  computer  program  which  evaluated  Welbull 's  formula  for 
fast  fracture  reliability,  which  In  turn  was  input  to  another  computer 
program  which  calculated  time  dependent  reliability.  From  double  tor¬ 
sion  testing  the  crack  velocity  exponent,  the  preraultipller  "A",  and 
the  critical  stress  Intensity  factor  were  obtained.  With  these  para¬ 
meters  and  the  value  of  the  str^^s  intensity  factor  coefficient,  "Y", 
the  value  of  "B"  was  calculated  .  The  value  of  "B"  and  crack  velo¬ 
city  exponent  “n"  were  input  to  a  computer  program  that  evaluated  the 
reliability  versus  time. 


FAST  FRACTUfiC  DOUBLE  TORSION 

testing  testing 


REUABILITT  VS  TIME 


Figure  15  -  Calculation  flow  path  utilizing  double  torsion  data 


Figure  16  shows  the  calculation  flow  path  utilizing  stress  rate 
testing  data  to  obtain  reliability  versus  time.  Stress  rate  testing 
determined  pairs  of  fast  fracture  strength  and  the  corresponding  value 
of  stress  rate.  The  value  of  "B"  and  the  crack  velocity  exponent  "n" 
were  calculate^g  pairs  of  strength  and  stress  rate  with  the 
formulas  shown  ,  .  The  fast  fracture  calculation  was  tlie  same  as 
that  used  with  double  torsion  data. 


FAST  FRACTURE  STRESS  RATE 

TESTING  TESTING 


RELIABILITY  VS  TIME 


Figure  16  -  Calculation  flow  path  utilizing  stress  rate  data 

The  third  data  base  available  for  predicting  time  decendent 
rellabi:^ity  Is  stress  rupture  testing.  Quinn^  ,  Trantina",  and 
Govlla  have  done  extensive  stress  rupture  testing  on  NC-132.  An 
empirical  formula  suitable  for  the  calculation  of  time  dependent 
reliability  using  stress  rupture  data  is  given  in  Appendix  A.  The 
calculation  flow  path  is  shown  in  Fig.  17.  The  reliability  versus 
time  formula  was  evaluated  for  each  element  in  the  finite  element 
mesh.  The  fast  fracture  reliability  calculation  is  the  same  as  that 
used  for  the  double  torsion  and  stress  rate  data  bases. 

The  procedures  outlined  here  were  used  to  design  the  disc.  The 
design  was  arrived  at  by  checking  trial  designs  with  the  reliability 
versus  time  calculation  procedures  until  a  satisfactory  design  was 
obtained.  The  reliability  versus  time  for  the  final  design  and  test 
conditions  is  shown  in  Fig.  18.  Note  the  widely  different  results  for 
the  different  data  bases. 
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RELIABILITY  (PERCENT) 


FAST  FRACTURE  STRESS  RUPTURE 
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RELIABILITY  VS  TIME 


Figure  17  -  Calculation  flow  path  utilizing  stress  rupture  data 


TIME  (HOURS) 


Figure  18  -  Calculated  reliability  versus  time. 
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III.  DISC  MACHINING 


Life  prediction  discs  were  machined  from  six  inch  square,  one  and 
three-eighths  inch  thick  billets  of  hot  pressed  silicon  nitride 
(NC132).*  The  four  inch  diameter  discs  were  made  in  several  steps 
beginning  with  a  slicing  operation  in  which  an  octagonal  block  was  cut 
out  of  each  billet.  This  piece  was  ground  to  the  approximate 
thickness  of  the  finished  disc  and  a  center  hole  or  bore  was  then  core 
drilled  through  with  a  diamond  core  drill.  The  bore  was  then  ground 
to  finish  size  allowing  the  disc  to  be  centrally  mounted  for  further 
grinding  operations.  Most  of  the  grinding  was  done  on  a  Brown  and 
Sharp  Universal  grinder.  Initial  material  removal  was  accomplished 
using  120-150  grit  diamond  plate  wheel  and  a  commercial  water  soluble 
coolant  (X50-2  Chem-trend).  Final  dimensions  were  produced  with 
200-250  grit  wheels.  The  disc  face  contours  were  produced  by  specially 
shaped  contour  grinding  wheels  as  shown  in  Fig.  19.  These  wheels  were 
rotated  at  12-1300  rpm  and  were  fed  into  the  disc  side  In  a  three  step 
operation  as  shown  in  Fig.  20.  The  disc  was  rotated  at  100  rpm  in  the 
opposite  direction  to  the  grinding  wheel  which  allowed  control  of  the 
grinding  and  stability  of  the  setup.  Finish  grinding  of  the  disk 
face  removed  .0001  inch  of  material  per  revolution  while  advancing 
radially  at  a  rate  of  .0003  to  .0005  per  revolution.  Grinding  wheel 
speed  was  7000  rpm. 

All  grinding  tools,  Including  contour  wheels,  were  diamond  coated 
and  final  polishing  of  the  disc  contours  was  done  with  diamond  pastes 
on  a  brass  lap.  The  lap  was  in  the  form  of  a  ball  on  a  1/4”  diameter 
shank  and  the  lap  was  driven  by  a  hand  held  electric  drill.  The  disc 
was  slowly  rotated  on  an  arbor  and  the  contour  was  polished  to  a  1-2  u 
inch  finish  by  successively  using  diamond  pastes  of  30-40,  6-12  and 
1-5  micron  grit. 

Finally  the  curvic  couplings  used  to  mount  the  disc  to  the  shaft 
were  ground  by  diamond  coated  grinding  wheels  on  special  purpose  cur¬ 
vic  grinding  equipment.  The  discs  were  then  Zygloed  for  possible 
cracks  caused  by  grinding  and  then  dimensionally  Inspected.  Figure  21 
shows  the  completed  disc  ready  for  testing. 


* 


Certified  Billets  from  Norton  Company 


Figure  19  -  Contour  Grinding  Wheel 
3  STEP  PROCEDURE  FOR  DISC  CONTOURING 


SURFACE  "A"  AND  RADIALLY 
INTO  BASIC  HUB  DIAMETER 


FINISHED 


STEP  3 _ 

GRINDING  WHEEL  ADVANCED 
PERPENDICULAR  TO  GRINDING 
WHEEL  AXIS  TO  COMPLETE 
CONTOURING. 


MATERIAL  IN  CROSS  HATCHED 
VOLUME  REMOVED. 


Figure  20  -  Procedure  for  disc  contouring 


Figure  21  -  Disc  ready  for  testing 

Figure  22  shows  a  lOX  photograph  of  the  disc  throat  region  ready  for 
testing. 


Figure  22  -  Disc  throat  region  at  lOX 
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rv.  TEST  PROGRAM 


All  Che  data  required  to  analyze  the  rig  performance  and  the  con¬ 
ditions  affecting  the  rotor  were  recorded  on  a  magnetic  tape  recorder. 
The  tape  had  a  four  hour  running  cycle  after  which  it  was  changed  and 
the  tape  erased  if  the  running  period  had  been  uneventful.  If  a 
failure  had  occurred,  the  magnetic  tape  was  played  back  into  a  Data 
Trak  recorder  on  to  light  sensitive  recording  paper.  The  resulting 
traces  gave  continuous  records  of  rig  speed,  disk  surface  temperature, 
cooling  air  flows  etc.  right  thru  a  failure.  The  channels  of  the  Data 
Trak  were  carefully  calibrated  before  testing  began  and  scaled  Mylar 
overlays  allowed  the  running  parameters  during  a  failure  to  be  very 
accurately  measured. 

These  records  invariably  showed  smooth  steady  state  running  con¬ 
ditions  right  to  the  very  Instant  of  disk  failure  at  which  time  an 
Instantaneous  change  in  disk  temperature  would  occur  followed  by 
slight  changes  in  oil  pressure,  cooling  air  flows  etc. 

The  disc  face  temperature,  displayed  on  the  tapes,  was  measured 
by  an  Ircon  radiation  pyrometer  which  viewed  the  face  of  the  spinning 
heated  rotor  thru  a  hole  in  the  face  plate  of  the  rig.  There  were  no 
intervening  glass  windows  or  other  obstructions.  This  unit  had  been 
calibrated  using  an  electric  furnace  containing  a  small  block  of  NC 
132  therraocoupled  with  Platlura  13%  Rhodium  couple  readout  on  a  Doric 
meter.  The  block  was  heated  to  a  temperature  higher  than  the  maximum 
of  the  pyrometer  range  being  calibrated,  the  furnace  heating  elements 
were  turned  off  and  readings  of  the  thermocouple  and . radiation  pyro¬ 
meter  were  taken  as  the  block  slowly  cooled  down  thru  the  range  being 
calibrated.  The  thermocouple  was  placed  inside  the  block  within  1/8” 
of  the  face  whose  temperature  was  being  measured  by  the  pyrometer. 
The  technique  of  calibrating  during  cooling  was  necessary  because  the 
thin,  white  hot  furnace  heating  elements  affected  the  pyrometer 
readings.  However,  when  the  power  was  shut  off,  the  heating  wires 
quickly  cooled  while  the  block  had  hardly  responded.  The  thermocouple 
had  a  factory  accuracy  of  ±  5®?  and  this  allowed  an  emlssivlty  setting 
of  .88  to  be  selected  which  gave  an  error  of  only  +10®F  at  the  test 
temperature  of  2300*F.  This  10®  correction  was  taken  into  account 
when  setting  the  temperature  of  the  disks  during  actual  testing  effec¬ 
tively  limiting  the  error  of  disk  temperature  measurement  to  -  5®F, 


Actual  testing  of  a  disk  was  begun  with  a  cold  spin  test  to  full 
test  speed  (50000  rpm)  to  assure  that  no  mechanical  problem  would 
require  a  shut  down  after  the  disc  had  been  brought  to  test 
temperature.  The  disc  was  then  spun  at  2000  rpm  and  the  gas  burners 
were  lit  and  the  disc  rim  allowed  to  heat  up  to  1800  -  2000®F  range. 
The  speed  would  then  be  Increased  to  10,000  rpm  and  the  rim  stabilized 
at  2100'2200®F.  A  second  speed  Increase  would  then  take  the  disc  to 
24000  rpm  and  a  15  minute  warm  up  allowed  during  which  bolt  cooling 
air  would  be  carefully  adjusted  and  disc  rim  would  be  brought  to 
2150-2200“F.  Holding  Che  rim  below  2300“F  before  reaching  50,000  rpm 
was  necessary  because  Increases  in  speed  had  the  effect  of  slightly 
raising  rim  temperature  all  other  conditions  remaining  the  same,  thus 
a  possible  over-tempera  e  condition  was  avoided. 

When  all  conditions  were  throughly  stabilized  at  acceptable 
levels,  the  speed  was  Increased  to  50,000  rpm  and  any  final  rim  tem¬ 
perature  adjustments  necessary  to  get  a  2300“F  temperature  were  made. 
When  the  rim  came  to  2300“F  at  50000  rpm  the  test  time  was  officially 
begun.  This  steady  state  condition  was  held  until  Che  disc  either 
failed  or  ran  for  24  hours.  A  disk  failure  was  a  barely  audible 
"thump"  after  which  the  rig  ran  smoothly  along. 

Post  failure  investigations  usually  showed  the  disc  as  well  as 
the  spacer  and  burst  ring  to  be  completely  shattered.  The  metal 
attaching  parts,  bolt  curvic  adaptor  and  nut,  were  usually  found  bro¬ 
ken  off  although  after  several  failures  the  bolt  was  found  in  place, 
straight  and  neither  stretched  or  over-heated.  This  was  not  unex¬ 
pected  since  the  bolt  had  originally  been  designed  for  a  stretch  load 
of  4700  lbs  which  took  it  to  80%  of  its  tensile  strength  at  running 
temperatures.  The  installation  load  of  only  2500  lbs  thus  left  a  very 
comfortable  safety  margin  and  the  critical  cooling  air  temperature  was 
measured  by  a  thermocouple  at  its  exit  from  the  bolt  making  precise 
control  very  easy. 


V.  EXPERIMENTAL  RESULTS 

A  total  of  twelve  discs  were  tested  in  the  program.  The  discs 
were  identified  by  the  four  digit  billet  serial  number  assigned  by  the 
Norton  Company.  Ford  Motor  added  the  prefix  NC  to  distinguish  these 
discs  from  others  fabricated  by  Ford  for  other  programs.  Two  discs 
failed  on  start-up  and  were  not  Included  as  part  of  the  database 
because  these  failures  were  probably  not  due  to  time  dependent  pheno¬ 
mena.  Disc  number  NC-1321  was  shut-down  after  thirteen  hours  of 
steady  state  testing,  but  it  failed  on  an  attempted  restart  so  it 
was  considered  censored  at  13  hours.  Disc  numbers  NC-1319,  and 
NC-1337  were  censored  at  24  hours  and  disc  NC-1323  was  censored  at  25 
hours.  These  discs  were  shut  down  after  continuous  steady  state 
testing  without  failure. 


The  test  results  are  summarized  in  Table  VI.  The  survival  proba¬ 
bility  of  each  disc  was  calculated  with  the  following  formula  given  by 
Bury^  . 


-  k 

R(X')  =  n 


i 


N-i 

N-i+1 


(4) 


The  format  of  the  table  is  also  taken  from  Bury.  The  experimental 
failures  are  plotted  on  Fig,  23  as  circles.  These  points  form  a  data 
base  against  which  analytical  methods  of  predicting  reliability  versus 
time  for  complex  structures  can  be  measured.  Also  shown  on  Fig.  23 
are  the  calculated  reliabilities  versus  time  from  the  three  data  bases 
used.  These  are  the  same  as  Fig.  18.  The  results  calculated  from 
stress  rupture  data  appear  to  give  the  best  correlation. 


F Lgure  23  -  Experimental  and  Calculated  Reliability  Comparison 


As  an  estimate  of  error  in  the  calculated  reliabilities  the 
radiation  pyrometer  was  assumed  to  have  an  error  of  -20°F  and  the 
results  were  recalculated  and  plotted  on  Fig.  24.  The  stress  rate 
error  was  too  small  to  appear  on  the  graph  but  the  calculated  reliabi¬ 
lity  at  25  hours  is  0.966,  .961,  and  .956  for  2280®F,  2300‘*F  and 
2320“F  maximum  material  temperatures  respectively.  The  actual  pyro¬ 
meter  temperature  error  was  estimated  to  be  ±5“F  because  it  was 
calibrated  with  a  sample  of  hot  pressed  silicon  nitride  instrumented 
with  a  thermocouple.  The  results  in  Fig.  24  serve  as  an  upper  limit  on 
the  experimental  and  analytical  determination  of  the  temperatures  and 
their  effect  on  ce Iculated  reliabilities. 


Figure  24  -  Effect  of  estimated  errors  in  temperature 

upon  calculated  results. 


VI.  CONCLUSIONS 

The  comparison  of  calculated  to  experimental  results  using  the 
three  different  databases  show  that  stress  rupture  data  is  the  best 
foL  predicting  reliability  versus  time. 


VII.  RECOMMENDATIONS 


Future  time  dependent  reliability  Investigations  should  con¬ 
centrate  on  stress  rupture  experimental  methods  and  analytical  tech¬ 
niques  to  predict  time  dependent  reliabilities  using  stress  rupture 
data. 


The  discs  from  the  test  program  were  delivered  to  AMMRC  for 
study,  and  it  is  recommended  the  AMMRC  conduct  a  f ractographlc 
investigation  of  the  discs. 
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APPENDIX  A 


Stress-Rupture  Formula 

Use  of  thCj^  t^wer  relation  of  crack  velocity  as  a  function  of 
stress  intensity  , 


ak: 


A-1 


12  14 

leads  tc  a  reliability  versus  time  equation  of  the  form  , 


2  -2- 

R  -  R^f  ''  B  '' 


A-2 


Examination  of  the  experimental  evidence  from  a  number  of  stress  rup¬ 
ture  testa  Indicate  that  the  results  fit  an  equation  of  the  form 


R  -  R^^  exp  (-<|i<t))  A-3 

where  (Kt)  is  some  empirical  relationship  of  the  variables  that 
describe  the  time  dependent  behavior  of  the  material.  For  a  constant 
temperature  these  variables  are  the  stress,  the  volume,  and  the  time. 
The  form  of  4>(t)  is  assumed  to  be 


<p(t) 


Veff„ 
_ 2  y 

Veff^ 


A-4 


The  subscript  1  applies  to  the  database  from  which  the  constants  C,  n, 
y,  oj,  and  tj  are  determined.  Subscript  2  applies  to  the  element  for 
which  reliability  is  being  calculated. 


( 


;  i 


'  1 


•  .t 


fa 


The  crack  velocity  exponent  was  chosen  because  stress  rupture  tests  in 
magnesia  doped  silicon  nitrides  are  described  by 

o"t  »  constant  A-5 


when  plotted  as  stress  versus  time.  The  effective  volume  was  used  as 
taken  from  Welbull  theory  to  correct  for  size  differences  between  the 
elements.  Equation  A-4  may  be  modified  to  use  the  characteristic 
strength  of  the  st^ycture  by  substituting  the  following  equation  given 
by  Paluszny  and  Wu  . 


V 


eff 


A-6 
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since  at  the  same  temperature  Cqj 
0  *^0  t„ 

Kt)  =  C(-i)"  (  — )^  " 

1  *^1 


<3q2  equation  A-4  becomes 

A-7 


For  a  MOR  bar  Oq  is  the  characteristic  modulus  of  rupture,  MOR,  which 
makes  the  application  of  Equation  A-7  relatively  simple.  Substitution 
of  Equation  A-7  into  Equation  A-3  yields  the  form  of  the  stress  rup¬ 
ture  equation  used. 


0  6  t 

R  -  R  ,  exp  <-C(^)"  )  A-8 
1  6^  1 

Equation  A-8  is  for  one  temperature.  Data  was  available  at  three  tem¬ 
peratures  1100°C,  1200®C,  and  1300*C.  For  intermediate  temperatures 
Oqj^  n,  and  m  were  linearly  Interpolated. 

To  apply  this  equation  to  the  disc,  the  characteristic  strength, 
Og-,  of  each  element  was  calculated,  as  well  as  the  fast  fracture 
reliability,  Rff.  The  reliability  for  the  entire  disc  was  calculated 
from 


R-R 
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R^ 


«  «  « 


R  . 

q 


A-9 


The  form  of  Equation  A-8  makes  it  less  sensitive  to  the  calcula¬ 
tion  of  fast  fracture  reliability  than  Equation  A-2.  In  equation  A-2 
the  fast  fracture  reliability  is  raised  to  a  power;  therefore,  the 
result  is  very  sensitive  to  the  fast  fracture  reliability.  Equation 
A-4  has  a  volume  correction  in  its  time  dependent  term  which  appears 
to  be  necessary  since  ceramics  are  proven  to  be  sensitive  to  volume  in 
fast  fracture  and  this  should  be  true  in  static  fatigue.  Temperature 
dependence  can  be  Included  in  equation  A-4  by  the  use  of  another  term 
so  that  4>(t)  become 


Oo  eff„  T-  t, 

<t>(t)  -  c  (- — ^ 

°1  eff^  n  ^2 


A-10 


This  was  not  done  in  this  study  because  of  the  lack  of  sufficient  data 
to  determine  z.  In  this  study  the  temperature  dependence  was  handled 
by  considering  n,  ,  and  y  as  functions  of  temperature. 
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To  utilize  the  empirical  stress  rupture  formula  given  here 
requires  a  body  of  stress  rupture  data  from  which  the  empirical 
constant|  can  be  determined.  The  data  base  was  Laken  from  Quinn  and 
Trantlna  and  the  results  are  given  in  Table  The  reference 
stress  was  taken  directly  from  Quinn's  curves  of  stress  versus  time 
using  time  equals  one  hour,  the  reference  time.  2?^®  stress  given  by 
Quinn  was  multiplied  by  ,82  to  account  for  creep  .  The  crack  velo¬ 
city  exponents  were  taken  from  Quinn's  graph.  The  reference  charac¬ 
teristic  MOR  is  the  expected  MOR  scaled  to  the  size  bar  used  by  Quinn. 
The  MOR  exponent  was  determined  by  the  amount  required  to  correct  for 
the  difference  in  bar  size  and  loading  configuration  between  Quinn's 
and  Trantina's  tests.  They  used  different  effective  volumes. 

The  formula  does  not  necessarily  give  good  predictions  in  all 
cases  as  illustrated  ii^^Fig.  A-1.  This  shows  a  prediction  for  stress 
rupture  data  of  Govila  taken  on  a  different  size  bar  than  Quinn  or 
Trantina  used.  The  result  is  poor.  One  possible  explanation  may  be 
billet  to  billet  variation  of  the  material.  The  applicability  of  the 
proposed  empirical  formula  requires  further  investigation;  however,  it 
was  proposed  in  this  report  as  the  method  that  gave  the  best  result. 


TIME  (HOURS) 


Figure  A-1  -  Experimental  and  Calculated  Reliability 
Comparison  For  415  MPa  Stress  Rupture  Tests, 
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APPENDIX  B 


Calculation  of  Constant  "B" 

The  calculation  of  time  dependent  reliabilities  is  based  on  the 
power  relation  of  crack  velocity  as  a  function  of  stress  intensity 
factor. 


V  =  akJ  B-1 

The  values  of  A  and  n  are  determined  experimentally,  so  a  literature 
search  was  conducted  to  find  published  values  to  utilize  for  design 
and  analysis.  Table  III  gives  the  experimental  values  that  were 
found.  The  experimental  values  were  published  in  graphical  or  numeri¬ 
cal  form  and  the  manner  of  presentation  Is  shown  on  the  table.  In  the 
cases  of  graphical  publication  the  numerical  values  were  found  by 
taking  values  of  velocity,  *'V”,  and  stress  intensity  factor,  "Kx", 
from  the  graph  and  substituting  into  Eqn.  B-1.  It  should  be  noted 
that  considerable  error  in  the  determination  of  A  and  n  may  result 
from  this  procedure.  Cargill  and  Annis  published  their  results 
graphically;  however,  they  furnished  the  values  of  A  and  n  determined 
by  least  sq^res  flttin^.^of  their  data.  The  A  and  n  values  of  Annis 
and  Cargill  and  Govlla  were  used  in  this  prograii. . 

The  values  of  A  and  n  available  are  published  for  930“C,  1200“C, 
and  1j00“C.  Reliability  calculation  requires  values  for  A  and  n  at 
the  centroidai  temperature  of  each  finite  element  of  the  finite  ele¬ 
ment  grid  shown  in  Fig,  2.  To  obtain  the  required  values,  linear 
interpolation  was  used. 

The  critical  stress  intensity  factor  for  the  material  is  an 
Important  parameter  when  double  ^grslon  data  is  used  to  calculate  time 
dependent  reliability.  Larsen  presented  graphical  data  for  hot 
pressed  silicon  nitride  and  the  numerical  values  used  are  shown  in 
Table  IV,  Also  required  for  the  use  of  double  torsion  data  is  an 
assumption  about  the  shape  of  the  flaw,  so  that  the  stress  intensity 
factor  coefficient  may  be  determined.  For  this  program  the  flaws  were 
assumed  to  be  penny  shaped  internal  flaw  with  a  stress  intensity  fac¬ 
tor  coefficient'^  of  2//Tr, 

As  shown  in  Fig.  15  the  value  of  ”B'’  must  be  computed  before 
reliability  versus  time  can  be  calculated  and  the  following  equations 
show  the  evaluation  of  "B"  at  1200®C  with  the  required  material  and 
geometric  constants  and  their  sources. 


References 


n 

«  27.4 

27 

A 

,  ,p,-23  ra 

=  1.58x10  - 

sec 

1 

(MPa/ra)^ 

27 

^IC 

*7.2  MPa/m 

28 

Y 

»  2//x 

29 

B 


_ 2 _ 

(27.4-2)(^)2  1.58  X  10"^^  x  (7.2)^^*^"^ 


hr. 

3600  sec 


B  -  0.000182  MPa^  -  hr. 

An  alternate  data  base  to  use  for  predicting  time  dependent 
failure 2^3  stress  rate  data.  The  NC-132  stress  rate  data  Is  given  In 
Table  V  .  Reference  26  gives  data  for  several  stress  rates,  but  the 
lowest  available  stress  rate  and  corresponding  fracture  strength  was 
chosen  on  the  assumption  that  slow  stress  rate  result  was  most  appli¬ 
cable  for  predicting  stress  rupture  results.  As  shown  on  Fig.  16  the 
value  of  "B"  must  be  evaluated  before  reliability  versus  time  can  be 
calculated  and  the  following  equations  show  a  typical  calculation. 
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Temperature  =  1204°C 

Reference 

n 

= 

19.9 

26 

a 

485  MPa 

26 

• 

0 

= 

0.8  MPa/min. 

26 

B 

a(n4-l ) 

B 

(485)^ 

0.8x( 19.9+1) 

hr. 

*  60  min. 

B 

-  113720  MPa^  • 

-  hr. 

APPENDIX  C  NOMENCLATURE 


A  Premulttplier  in  crack  velocity  equation,  meter/second 
(MPa-meter)'^. 

B  Constant  in  reliability  versus  time  equation,  MPa'’-hr. 

C  Constant  in  stress  rupture  equation, 

1  Order  of  a  specimen, 

k  Index  on  Operator  II. 

Kj.^  Critical  stress  intensity,  a  material  parameter,  MPa/meter 
Kj  Stress  intensity  Mpa/meter. 

m  Weibull  modulus. 

MORq  Characteristic  modulus  of  rupture,  MPa. 
n  Crack  velocity  exponent. 

N  Initial  sample  size. 

(t»(t)  Function  in  stress  rupture  equation, 

q  Number  of  elements  in  the  structure, 

R  Reliability  as  a  function  of  time. 

Tff  Fast  fracture  reliability. 

R(X’)  Maximum  likelihood  estimator  of  reliability. 

0  Stress,  MPa. 

Og  Characteristic  strength  of  the  structure,  MPa. 

Fracture  stress,  MPa. 

a  Stress  rate,  MPa/min. 

3 

Og  Weibull  parameter,  MPa/ (meter  )  l/m, 

T  Temperature,  ®C. 

t  Time,  hours. 

3 

V  Volume,  meter  . 

V  Crack  velocity,  meter/sec. 

3 

Vgff  Effective  volume,  meter 

V  Stress  intensity  factor  coefficient,  non-dimensional, 

y  Exponent  on  effective  volume. 

z  Exponent  on  temperature. 


-32- 


TABLE  T 


Thermal  Properties 


Temperature 

<‘F 

Thermal 

Conductivity 

Btu 

hr  ft®F 

Specific 

Heat 

Btu/lb„-®F 

Spin  Disc  and  Ceramic 

70. 

17 

.18 

Curvic  Spacer  Material 

500. 

15. 

.23 

Hot  Pressed  Silicon  Nitride  1000. 

13. 

.26 

Density  =■  3.18  gm/cra 

1500. 

11. 

.29 

2000. 

9.2 

.33 

2500. 

8.0 

.32 

Nut  Material 

0. 

6.0 

.10 

Inconel  X-750  ^ 

Density  =■  8.30  gm/cm 

3000, 

20.5 

.18 

Tie  Bolt  Material 

0. 

6.0 

.10 

Inconel  X-718  ^ 

Density  ■  8.19  gm/cm^ 

3000. 

21.8 

.10 

Compression  Member 

0. 

6.5 

.11 

Material 

A-286  ^ 

Density  ■7.91  gm/cm’ 

3000. 

25.2 

.11 

Shaft  Material 

0. 

21.7 

.11 

AMS  6265  3 

Density  »  7.83  gm/cm’ 

3000. 

21.7 

.11 

Curvic  Material 

70. 

3.99 

.18 

500. 

3.87 

.23 

Density  *  3.18  gm/cm 

1000. 

3.72 

.26 

1500. 

3.54 

.29 

2000. 

3.33 

.33 

2500. 

3.16 

.32 

Temperature 

Thermal 

^ecif  ic 

Density 

op 

Conductivity 

Heat 

Btu 

Btu 

Ihig 

hr  ft  “F 

Ibrn^F 

ft* 

Air  0 

,013 

.239 

.0863 

100 

.015 

.240 

.0680 

500 

.023 

.247 

.0392 

800 

.028 

.256 

.0299 

1000 

.032 

.262 

.0258 

1500 

.040 

.27  6 

.0192 

2000 

.047 

.286 

.0153 

2500 

.051 

.292 

.0127 

3000 

.054 

.297 

.0108 
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TABLE  II 


Elastic  Properties 


NC-132 

Hot  Pressed 
Silicon 

Temperature 

Young 's 
Modulus 

Poisson's 

Ratio 

Shear 

Modulus 

Coefficient 
of  Thermal 
Expansion 
ln/in/°F 

Density  ■  3 

“F 

0 

xlO  psi 

0 

xlO  psi 

xl0"6 

3.18  gm/cm 

78 

41.9 

.219 

17.’ 

.7  4 

500 

41.9 

.203 

17.5 

.97 

1000 

41.8 

.190 

17.5 

1.32 

1500 

40.6 

.189 

17.1 

1,53 

2000 

38.7 

.194 

16.1 

1.68 

2500 

35.4 

.194 

14.3 

1.82 

Reference  13 


TABLE  III 


Crack  Velocity  Parameters  Survey 


Material 


Crack 
Velocity 
Exponent 
n 


HPSN 

NC-132 

930 

6.18x10-56 

HPSN 

NC-132 

1200 

1.58x10-23 

HPSN 

NC-132 

1300 

4.3x10-9 

HPSN 

NC-132 

1350 

7.9x10-9 

HPSN 

NC-132 

1400 

1x10-7 

HPa^ 

NC-132 

1200 

1.21x10-31 

HPSN 

NC-132 

1400 

9x10-8 

HPSN 

HS-130 

1250 

1.41x10-1^ 

HPSN 

HS-130 

1300 

9.31x10-1^ 

HPSN 

HS-130 

1350 

9.83x10-13 

HPSN 

HS-130 

1400 

3.05x10-12 

HPSN 

HS-130 

1200 

2.56x10-59 

HP  a? 

HS-130 

1400 

8.22x10-8 

Form  of 
Data 

Presentation 


Graphical 

Graphical 

Numerical 

Numerical 

Numerical 

Graphical 

Graphical 

Graphical 

Graphical 

Graphical 

Graphical 

Graphical 

Graphical 

Numerical 


TABLE  IV 


Critical  Stress  Intensity  Factor 


Temperature 

®C 


18 

930 

1200 

1300 

Reference  28 


Stress  Intensity 
Factor 
MPa/m 


5. 

5, 

7.2 

11.0 


TABLE  V 


Stress  Rate  Data  of  Norton  NC-132  HPSN 


Temperature 

"C 

Characteristic 
MOR 
’  MPa 

Welbull 

Modulus 

Stress 

Rate 

MPa /min. 

Crack 

Velocity 

Exponent 

Number 

of 

Samples 

21 

787. 

6.5 

1870.0 

72.2 

30 

704 

730. 

6.0 

1.8 

72.2 

5 

871 

818. 

8.5 

1.8 

72.2 

5 

1038 

666. 

9.9 

1.5 

72.7 

5 

1204 

485. 

8.0 

.8 

19.9 

5 

1371 

313. 

11.2 

91,0 

17.2 

15 

Reference  26 
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TABLE  VI 


Summary  of  Discs  Tested 


Discs 

Number 

Order 

Number 

1 

Failure  Censoring 
Time  Time 

Hours  Hours 

N-1 

N-i+1 

Reliability 

ft(Xi) 

NC-1333 

1 

.20 

9/10 

.90 

NC-1322 

2 

7.17 

8/9 

.80 

NC-1318 

3 

8.50 

7/8 

.70 

NC-1321 

4 

13.00 

- 

NC-131A 

5 

13.83 

5/6 

.58 

NC-1335 

6 

18.27 

4/5 

.47 

NC-1324 

7 

18.58 

3/4 

.35 

NC-1319 

8 

24.00 

NC-1337 

9 

24.00 

NC-1323 

10 

25.00 

TABLE  Vir 


Stress  Rupture  Data  Base 


Temp. 

Constant 

Reference 

Stress 

MPa 

Crack 

Vel. 

Exp. 

Reference 

Char. 

Mpa 

MOR 

Exp. 

Reference 

Time 

Hours 

1100 

.69 

402 

19 

634.8 

8.15 

1. 

1200 

.69 

295 

12 

517.0 

4.80 

1. 

1300 

.69 

230 

10 

411.1 

4.50 

1. 

Reference  24 
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20000  Rotunda  Drive,  Dearborn,  MI  48121 
1  ATTN:  Mr.  A.  F.  McLean 

1  Mr.  £,  A.  Fisher 

1  Mr.  J.  A.  Mangels 

1  Hr.  R.  Govila 

General  Electric  Company,  Research  and  Development  Center, 

Box  8,  Schenectady,  NY  12345 
1  ATTN:  Dr.  R.  J.  Charles 

1  Dr.  C.  D.  Greskovich 

1  Dr.  5.  Prcchazka 

General  Motors  Corporation,  AC  Spark.  Plug  Division, 

Flint,  Ml  48556 
1  ATTN:  Or.  M.  Berg 

Georgia  Institute  of  Technology,  EES,  Atlanta,  GA  30332 
1  ATTN:  Mr.  J.  D.  Walton 

GTE  Lab.nratories,  Waltham  Research  Center,  40  Sylvan  Road,  Waltham, 
MA  02154 

1  ATTN:  Dr.  C.  Quackenbush 

1  Dr.  w.  H,  Rhodes 

I  IT  Research  Institute,  10  West  35th  Street,  Chicago,  II.  60616 
1  ATTN:  Hr.  S.  Bortz,  Director,  Ceramics  Research 

Institul  fur  Werkstoff-Forschung,  OFVLR.  505  Porz-Wahn,  Linder 
Hohe,  Germany 
1  ATTN:  Or.  W.  Bunk 

International  Harvester,  Solar  Division,  2200  Pacific  Highway, 

P.O.  Box  80966,  San  Diego,  CA  92138 
1  ATTN:  Or,  A.  Metcalfe 

1  Ms.  M.  E.  Gulden 


Kawecki  Berylco  Industries,  Inc.,  P.O.  Box  1462,  Reading, 

PA  19603 

1  ATTN:  Hr.  R.  J,  Longenecker 

Martin  Marietta  Laboratories,  1450  South  Rolling  Road,  Siltiinoro, 

HD  21227 

1  ATTN:  Dr,  J.  Venables 

Massachusetts  Institute  of  Technology,  nepartment  of  Metallurgy 
and  Materials  Science,  Cambridge,  MA  02139 
1  ATTN;  Prof.  R.  L.  Coble 

1  Prof.  H.  K.  8owen 

1  Prof.  W.  0.  Kingery 

Midwest  Research  Institute,  425  Volker  Boulevard,  Kansas  City, 

HO  64110 

1  ATTN:  Mr.  Gordon  w.  Gross,  Head,  Physics  Station 

Norton  Company,  Worcester,  MA  01606 
1  ATTN:  Or.  N.  Ault 

1  Dr.  M.  L.  Tort! 

« 

Pennsylvania  State  University,  Materials  Research  Laboratory, 

Materials  Science  Department,  University  Park,  PA  16802 
1  ATTN:  Prof.  R.  E.  Tressler 

1  Prof.  R.  Bradt 

1  Prof.  V.  S.  Stubican 

R1A$,  Division  of  the  Martin  Company,  Baltimore,  MD  21203 
1  ATTN:  0»*.  A.  R.  C.  Westwood 

Stanford  Research  International,  333  Ravenswood  Avenue, 

Menlo  Park,  CA  94025 
1  ATTN:  Or.  P.  Jorgensen 

I  Dr.  D.  Rowel iffe 

State  University  of  New  York  at  Stony  Brook,  Department  of  Materials 
Science,  Long  Island,  NY  11790 
1  ATTN:  Prof.  Franklin  F.  Y.  Wang 

United  Technologies  Research  Center,  East  Hartford,  CT  06108 
I  ATTN:  Or.  J.  Brennan 

1  Dr.  F,  Galdsso 

University  of  California,  Lawrence  Livermore  Laboratory, 

P.O.  Box  808,  Livermore,  CA  94550 
1  ATTN:  Dr.  C.  F.  Cline 

University  of  Florida,  Department  of  Materials  Science  and  Engineering, 
Gainesville,  FL  32601 
1  ATTN:  Or.  L.  Hench 

University  of  Newcastle  Upon  Tyne,  Department  of  Metallurgy  and 
Engineering  Materials,  Newcastle  Upon  Tyne,  NEI  7  Ru,  England 
I  ATTN;  Prof.  K.  H.  Jack 

University  of  Washington,  Ceramic  Engineering  Division,  FB-lO, 

Seattle,  WA  98195 
1  ATTN:  Prof.  James  I.  Mueller 

westinghouse  Electric  Corporation,  Research  Laboratories, 

Pittsburgh,  PA  15235 

1  ATTN:  Or.  R.  J.  Bratton 

Director,  Army  Materials  and  Mechanics  Research  Center, 

Watertown,  MA  02172 

2  ATTN:  ORXMR-PL 

1  ORXMR-PR 

1  ORXHR-K 

10  ORXMR-MC,  Or.  R.  Katz 

10  ORXMft-SM,  Dr.  £.  Lenoe 


